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Abstract 

Based on an extended NJL model including heavy quark flavors, we calculate the form factors 
of pseudo-scalar and vector mesons. After taking into account the vector-meson-dominance effect, 
which introduces a form factor correction to the quark vector coupling vertices, the form factors 
and electric radii of vr"’' and pseudo-scalar mesons in the light flavor sector fit the experimental 
data well. The magnetic moments of the light vector mesons p"*" and K*~^ are comparable with 
other theoretical calculations. The form factors in the light-heavy flavor sector are presented to 
compare with future experiments or other theoretical calculations. 


PACS numbers: 12.39.Fe, 12.39.Hg, 14.40.-n 

Keywords: NJL model, heavy meson, form factor, magnet moment 


* dwz@pku.edu.cn 


1 





I. INTRODUCTION 


,Q 


The Nambu-Jona-Lasinio (NJL) model [ll, has been widely used in hadron physics as 
an effective model to study chiral symmetry in the degree of quark freedom. Usually this 
model deals with light hadrons composed of only light quark flavors u, d, s with SUf{3) 
symmetry 

In hadron systems including heavy flavors, such as light-heavy mesons, although the 
chiral symmetry is broken due to the mass of the heavy quark, a complementary heavy flavor 
symmetry emerges and the so-called heavy quark effective theory (HQET) was formulated 


141. In Ref. 151, the NJL model was 


for this using the technique of I/uiq expansion [?- 
extended to include heavy quark flavors to investigate such light-heavy mesons like D^*'> and 
RW mesons. 


In our previous work [^, we also tried to extend the NJL model to include heavy flavors 
by expanding the NJL interaction strengths in the inverse power of constituent quark masses 
according to HQET. Based on this extension, we obtained the meson masses and meson- 
quark coupling constants of all light and light-heavy mesons in a unified way. Furthermore, 
the decay widths of the mesons were calculated from those effective meson quark couplings 


izl. 


In this work, we further calculate the electromagnetic form factors of mesons within this 
extended model. Electromagnetic form factors play an important role in our understanding 


of hadronic structure, 
in several experiments 


he form factors of the pseudo-scalar mesons vr and K were measured 


18 


20| and in some previous theoretical works, the form factors of vr 


and K mesons were studied in the NJL model |2ll. |22||. After considering the effect of vector- 
meson-dominance of the vector mesons, such as the p meson, in the calculation, typically the 
form factors of vr fit the experimental data well. Furthermore, the form factor of tt was also 
studied in case of hnite temperature with the NJL model j^|. Certainly the form factor of 


vr was studied in many other theoretical approaches, such as the Dyson-Schwinger equa 


and the lattice QCD method 


using a conhning quark propa gator 1241 , light-cone or covariant quark wave functions 


25 


ion 




27 . 281. Also, with the QCD factorization approach [? ? ? ? 


], the form factor can be extrapolated to higher energy regions by taking into account the 
perturbative QCD contribution. 

The form factors of vector mesons have a rather more complicated structure. Con- 
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sequently they can provide us more information about vector mesons, such as magnetic 
moments and quadrupole moments. Presently there are only theoretical results about the 
form factors of vector mesons. Some works have used the constituent quark model and the 


light front dynamics 291431| or Dyson-Schwinger equations ^ . Lattice QCD calculation 

'unctions method [33[, and the background held 
35| . The magnetic moments of vector mesons 


have been performed with the three-point 
method using only two-point functions 34 


were also calculated by dynamics with the external magnetic held 
rules 

There are a few papers studying the form factor of light-heavy mesons 


36[ |. and with QCD sum 


These focus 


on the electroweak form factors. From the heavy havor symmetry, those form factor should 
be unifying described by the Isgur-Wise function when the heavy havor mass goes to inhnity. 


Here, we perform a systematic calculation of the meson form factors, including pseudo¬ 
scalar mesons and vector mesons, of both the light havor sector and the light-heavy havor 
sector, within the extended NJL model. In the next section, we will introduce our model 
and formalism. The numerical results and discussion will be presented in Section 3. 


II. MODEL AND FORMALISM 

A. Extended NJL model 

To deal with both light and heavy mesons in the Nambu-Jona-Lasinio (NJL) model, in 
Ref. the four-fermion point interactions are modihed to 

= GvmriWW) + + (?7„75A;<;)(9'7,75A;9')1 (1) 

nijqTTXqf 

where A“ are the generator of SU (3) in color space and q, q' = u, d, s, c, b including both the 
light and the heavy havors. Here the second part of the interaction is required to improve 
the spectra of light vector mesons and the factor of l/{mqmqi) guarantees that the symmetry 
of heavy havors will still be hold in the heavy quark limit according to HQET. 

By solving the Bethe-Salpeter equation (BSE), we obtain the meson masses and their 
coupling constants with quarks. We will use the ehective Lagrangian to describe the quark 
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FIG. 1. Feynman diagrams of meson form factor. 


interaction in mesons. In the case of vr and p, the effective Lagrangian reads 



( 2 ) 


Cpqq = - gpqQlpTq ■ p. 


( 3 ) 


Here the couplings g-^q^ g-^q and gpq are treated as constants since the energy of immediate 
quarks is truncated to the low energy region in the NJL model. 

In ref. [^, we have calculated the strong and radiative decays of vector mesons. In this 
work, we will use the above effective meson Lagrangian to further calculate the form factors 
of mesons. 


B. Form factor of pseudo-scalar mesons 

The dehnition of the form factor of a pseudo-meson is given by 


(7r+(p2)|V'7M^k+(pi)) = iPi+P2)pF{q‘^), 


( 4 ) 


where q = pi — p 2 is the transfer momentum. Its Feynman diagrams are shown in Fig. 1 
where mi and m 2 are the masses of the constituent quarks in the pseudo-scalar meson. 
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Using the Feynman rules, the amplitude reads 


{.Pl+P2)^,F{q^) ={Pi+P2)AQiP^^\q^) + Q2F^‘^\q^)] 

d^k 


(Pi +P 2 )^,F^^\q^) = - Tr 


(2vr) 


+ pi)i{g - g 


mi + m 2 


)i'lbS2{k) 


X ^9 + 9 -y- )ilbSi{k + P 2 ), 

mi + m 2 


( 5 ) 


( 6 ) 


(Pi = - Tr 


d'^k 


-f^,S2{k - p2)i{g + g 


X i{9-9- 


mi + m 2 
)*755'2(fc -Pi), 


)i'y5Si{k) 


(7) 


mi + m 2 

where and are the form factors of the quark and anti-quark respectively, Qi is the 
electron charge of i-th quark, 

S^{p) = - -^ (8) 

- m* + *e 

is the propagator of the Fth quark, and g and g are the coupling constants of the pseudo¬ 
scalar meson obtained in our previous work Q. 

In the Breit frame, P 2 ~ Pi = 0 pi = —p 2 - We introduce 


Pi =P + 9/‘^, 


P 2 =p-q/2, 


(9) 


where p = |(pi+P 2 ) = (Pi, 0), g = (0, q). Taking the direction of the z-axis along momentum 
Pi, we hnd 


(1) . f d^k \ 2^1 gg S2 gg S3 , g^ 

[q ) tUcUf j ^ mi + m 2 F mi + m 2 F ^ (mi + ms)^ F \ 

F*'^^(g^) =F*^^^ ^mi y-)- m2, q —)■ —q,p —pj, 


( 10 ) 

( 11 ) 


where 


S'! =S'o(mi,mi,m2), (12) 

5 2 =miS'o(mi, {k + p- g/2)^/mi, m 2 ) m25'o(mi, mi, /c^/m2), (13) 

5 3 =miS'o((fc + P + g/2)^/mi, mi, m 2 ) m25'o(mi, mi, k‘^/m 2 ), (14) 

5 4 =mlSo{{k + P + qjFf jmi, {k ^ p - g/2)^/mi, m 2 ) /c^S'o(mi, mi, m 2 ) 

-y mim2S'o((fc + P + g/2)^/mi,mi, k‘^/m2) 

+ mim2S'o(mi, {k + p - g/2)^/mi, k‘^/m 2 ) 

D =[(k + p + g/2)^ — mi -I- ie)[{k + p — g/2)^ — m\ + ie](k‘^ — m\ + ie). 
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(15) 

(16) 















and 


S'o(mi, m 2 , m 3 ) = 2 [m 3 (mi + m 2 ) - 2k ■ {k + p)] 


^ ‘ P r / \ ^2 1 ‘1 

+ 2—— m3(mi + m 2 ) - mim2 + p -k 

p2 


iVi]- 


(17) 


Note that the denominator D of the integrand is invariant nnder the transformation k —)■ —k. 
The electromagnetic radins will be farther obtained from the derivative of the form factor 


via 


We have 


where 


is the radins of the i-th. qnark. 



(18) 


(19) 

( 20 ) 


C. Form factor of vector mesons 


The dehnition of the form factor of the vector meson reads 


39 


40| 


(p (P2, A2)|V'7/2'0 Ip (Pi, Ai)) = - e*(p2, A 2 ) • e(pi, Ai)(pi +p2)M^i(g ) 

+ ^ 1 )? • ^*(P2, A2) - e* (p2, Aa)? • e(pi, \i)]F2{q^) 


q • e*{p 2 ,X 2 )q • e(pi, Ai) 

2m? 


(pi + P2)fF3{q^), 


( 21 ) 


where e(pi) and e{p 2 ) are the polarization vectors of the initial and the final vector meson 
respectively. Based on the Feynman diagrams, the LHS of Eq. fl^T]) can be written as 


ei/(pi, Ai)e^(p2, Aa)^))^, (22) 

where 

Gf =OiGO>“" + Q^Gf'\ (23) 

/ d^k 

+ Pi)'tgvrS2{k)igvpSi(k+p2), (24) 

- P2)igvl^Si{k)igv7'"S2{k - pi). (25) 
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Still working in the Breit frame and taking the z-axis along the momentnm pi, the 
polarization vectors are chosen to be 


e(Pi,±) =;^(0, l,±f,0). 

e(Pi,0) = —(pn,0,0,pio), 
m 


e(P2,±) =^(0, l,Ti,0,) 

e(P2,0) = —(p22,0,0,p2o). 
m 

(26) 


To retrieve Fi, we take the time component in Eq. and hnd that 


Ai)e^(p2, A2 )Go^ = - e*(p2, A 2 ) • e(pi, Ai)(pi + p2)oi^i(g^) 

g • e*(p 2 , Aa)? • e(pi, Ai) 


2m2 


-(Pi +P2)o^3(g )• 


Then Fi can be obtained via the transverse polarization 


(27) 


e^(pi,±)eA(P2,±)G;^ = -e'(P2,±) ■ e(pi,±)(pi + P2)on(?^). (28) 

To retrieve Fa, we take the spatial components in eq. fl^T]) and hnd that 

e;.(pi, Ai)e*(p 2 , A 2 )Gf'' = - [ei(pi, Ai)q ■ e*(p 2 , Aa) - e-(pa, Aa)^ ■ e(pi, Ai)]F 2 (g^) 

= {[e(pi, Ai) X e*(p 2 , Aa)] x qr}. F 2 (g^). (29) 

Still each form factor Fj is a charge weight average of form factors of qnark and anti-qnark 
in the vector meson, 

F,(g2) = Q,Fl^\q^) + Q2F\^\q^), (30) 

and 

= Pf' •H- ma, g -)■ -q,p -pj • 

Explicitly we obtain 

=«f(Pi, +)el(P2, -)Gj‘’‘'p(2po) 

2 f d‘k Gi 

-m,n,gvj (31) 

where 


Gi =4:[{k + p) ■ k - 17117712 + kl + ky] 

- 2^-^[p^ - k"^ - g^/4 -m\ + 2mim2 - 2{kl + kl)], (32) 

pZ y 
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FIG. 2. Feynmann diagram of the loop correction 
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2 


(33) 


where 


G 2 =4[fc ■ {k + p) - 77111712 + kl] 

~ ~ + 2(fc^ + kl)]. (34) 

We will not consider the form factor in this work. 


D. Vector-Meson-Dominance and Qnark Loop Correction 

According to the vector-meson-dominance pictnre, the tt and K form factor are dominated 
by the p, oj and 0 intermediate vector meson states [2l|. In the NJL model, the vector-meson- 
dominance is represented by the correction to quark-photon vertex as shown in the Feynman 
diagram in Fig. 2. The correction will introduce a form factor to the constituent quark j^. 
For the Ath quark 


F«(g2) = 




(35) 


where is the NJL vector coupling constant and J^y represents the transverse vector 
loop integral j4, 16|. The meson form factor will be modified to 


F{q^) = QiFW(g2)FW(g2) + 


(36) 
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TABLE I. The masses and quark coupling constants of pseudo-scalar mesons 



TT 

K 

D 

Ds 

B 

mass (Me V) 

139 

496 

1870 

1940 

5280 

9 

4.25 

4.32 

4.71 

5.03 

5.92 

9 

1.56 

1.61 

2.04 

2.09 

2.84 



TABLE 11. The masses 

and quark coupling constants of vector mesons 



P 

K* 

D* 

D*s 

B* 

mass (Me V) 

771 

918 

1990 

2120 

5310 

9 

1.29 

1.31 

1.64 

1.83 

2.51 


III. NUMERICAL RESULTS 


The parameters of the extended NJL model were fixed by fitting the meson mass spectra 
and decay constants in a previous work [^. The input parameters were the current masses 
of light quarks and the constituent masses of heavy quarks, two coupling constants and the 
3-dimensional cutoff: 


= 2.79MeV, 
•rric = 1.62GeV, 

A = O.SGeV, 
h = 0.65. 


= 72.0MeV, 
mb = 4.94GeV, 
Gv = 2.41, 


(37) 


Due to charge conservation, the form factor should be normalized to F{q‘^ = 0) = 1 for 
any hadron carrying -|-1 charge. We will make a self-consistent calculation by using the 
theoretical values of meson masses and the quark coupling constants together. This will 
guarantee the strict normalization of the form factor at = 0 23|. The theoretical values 
of pseudo-scalar and vector mesons are listed in Table [H and Table HIl respectively. 
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FIG. 3. The form factor of Tr"*" compared to the experimental data from Ref. 






FIG. 4. The form factor of K~^ compared to the experimental data from Ref. 



A. Pseudo-scalar mesons 

The form factors of tt’*' and are compared with experimental data in Fig. |3] and 
Fig. 0] respectively. The theoretical resnlts fit the experimental data well. In the theoretical 
calculation, the quark loop correction is included to account for the important effect of 
vector-meson-dominance. 

The heavy-light pseudo-scalar mesons like ^ and have no experimental data for 
form factor yet. In Fig. 01 we present the form factors of all positive pseudo-scalar mesons. 
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FIG. 5. The form factors of all positive pseudo-scalar mesons. 


TABLE III. Comparison of electric radii of pseudo-scalar mesons vs. experimental data. ri 

and r 2 are constituent quark and anti-quark radii in the meson. 


Meson 

(rf) (fm^) 

(r|) (fm^) 

(fm) 

(^^)exp. (fm) [12,^ 

TT^ 

0.322 

0.322 

0.57 

0.66 

K+ 

0.341 

0.206 

0.54 

0.56 

D+ 

0.080 

0.473 

0.46 


Dt 

0.083 

0.283 

0.39 


B+ 

0.795 

0.041 

0.74 



The form factor at low momentum can be well illustrated by the electromagnetic 
radius. The radii of all positive pseudo-scalar mesons are listed in Table IIIII 

As had been observed in Ref. 2lJ, it is the quark loop correction of vector-meson- 


dominance that makes the tt radius bigger than that of the K. In Eq. fll9p . the radius 
of a meson is a charge weight average of individual quark radii. From Eq. fl5^ . after con¬ 
sidering the quark loop correction, we have 


(r-) = (r^int + {r^)a , 


(38) 
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FIG. 6. The form factors of constituent quark and anti-quark in tt and K. 


where 




dq^ 


q2=0 


.dF, 


(i) 


dq^ 


. q2=o 


(39) 

(40) 


are the “intrinsic” charge radius and the quark loop correction respectively. The quark loop 
correction decreases as the quark mass increases, so the lighter quark has a larger radius 
than its heavier partner in any meson. We show the individual form factors of quark and 
anti-quark in tt and K mesons in Fig. [Bland also list the individual quark radii in Table HTTl 

Just like the and , the radii of the light-heavy mesons decrease 

as the meson mass increases. However the radius of the meson increases by roughly a 
factor 2. This mainly because, in a light-heavy meson, the heavy quark’s contribution is 
much smaller than that of the light one. If we neglect the contribution of the heavy quark, 
in the meson, the u-quark has a 2/3 charge weight of contribution. The d-quark, on the 
other hand, has only a 1/3 charge weight in the and mesons. The form factors of 
individual constituent quarks in D and B are shown in Fig. [71 
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FIG. 7. The form factors of constituent quark and anti-quark in B and D. 



B. Vector mesons 


The electric form factors Fi of vector mesons are shown in Fig. [HI The electric radii are 
listed in Table IIVI Becanse all vector meson masses are close to their thresholds, their bonnd 
energies are small and their radii are larger than their psendo-scalar partners. 

The magnetic form factors F 2 are presented in Fig. 0 They are connected to the magnetic 
momentnm through 


34| 


hv = 


^2(0)e 

2mv 


(41) 


The magnetic moments are also listed in Table IIVI The magnetic moments are given in the 
unit of nuclear magneton /i„. Generally, the magnetic momentum decreases as the meson 
mass increases. In our results, the magnetic moments of D* and are smaller than those 
of the light mesons p and K. However, the magnetic moment of B* is larger than that 
of D* and H*. The reason is still that the main contribution comes from the light quark 
but the charge of the u is larger than that of the d and s by a factor of 2. Up to now, no 


experimental data is avai 


other theoretical work 


able. We compare our results for the p"*" and K*^ mesons with 


29| and |33|. 
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FIG. 8. Fi of all positive vector mesons with respect to 


TABLE IV. Radii and magnetic moments of vector mesons 


meson 


rl{frri^) 

rdfm) 

Pl(Pn) 

h2(Pn) 

h(hn) 

h(hn)[29] 

p(Pn)[3^ 


1.267 

1.267 

1.12 

1.69 

0.85 

2.54 

2.56 

3.25 

K*+ 

1.304 

0.697 

1.05 

1.63 

0.63 

2.26 


2.81 

D*+ 

0.095 

1.366 

0.72 

0.42 

0.74 

1.16 



D*+ 

0.083 

0.567 

0.49 

0.42 

0.56 

0.98 



B*+ 

1.359 

0.025 

0.96 

1.4 

0.07 

1.47 




IV. SUMMARY 


With the extended NJL model including heavy flavors, we have made a systematic cal¬ 
culation of the form factors of mesons, including the pseudo-scalar mesons and their vector 
partners, of both the light flavor sector and the light-heavy flavor sector. The form factors 
of the 7T and K mesons £t the experimental data. Other form factors of mesons, especially 
of the light-heavy mesons, are presented here to compare with future experiments and other 
theoretical calculations such as lattice calculation. 
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